Tailoring the Fermi level of the leads in molecular-electronic devices 
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The dependence of the transport properties on the specific location of the Fermi level in molecular electronics 
devices is studied by using electrodes of different materials. The zero-bias transport properties are shown to 
depend dramatically on the elemental composition of the electrodes, even though the shape of the transmission 
coefficients is very similar. By using alkaline materials it is possible to move the Fermi level from the HOMO- 
LUMO gap to the LUMO resonance and change dramatically the length dependence of the conductance of 
molecular wires, which opens the possibility of using molecules with different lengths and very similar con- 
ductances in nanoscale circuits. This method shows how to dramatically increase the conductance of molecular 
devices and alter qualitatively and quantitatively their electronic and transport properties. 

PACS numbers: 85.65.+h,73.63.-b.73.40.-c 



INTRODUCTION 



During recent years, the race to develop a viable technol- 
ogy for sub- 10 nm electronic wires led a number of groups 
to study electron transport through single molecules bridg- 
ing two metalUc electrodes. Examples include both experi- 
mental (H H 1 H and theoretical S 0, S H works. In 
all of these experiments one or more of the metallic elec- 
trodes are formed either from break junctions ||2l |4|] or from 
the tip of an STM 111 |3|]. None of these contacting meth- 
ods is capable of being scaled to billions of devices on a sin- 
gle chip. One route to developing a scalable technology is 
to first deposit electrodes with a pre-defined gap and then to 
deposit molecules whose length is matched to the electrode 
gap. Since reproducible electrode gaps significantly below 
10 nm do not cuiTently exist whereas most studies of sin- 
gle molecules involve molecules of length < 3 nm, there is 
a need for systematic studies of families of molecules of vary- 
ing lengths which bridge those length scales. Current ex- 
amples include alkanedithiols and aUcanemonothiols lilQl . di- 
aminoacenes 11 1111 and oUgothiophenes lll2ll . 



When a low-enough voltage is applied to a single molecule 
bridging two metallic electrodes, the measured electrical con- 
ductance G is governed by the spatial extent of the molecular 
orbitals, their level broadening F due to the contacts and the 
position of the HOMO (-Eh) and LUMO (&0 levels relative 
to the Fermi energy (Ep) of the electrodes lll3[ll4ll5n . This 
dependence is captured by the Landauer formula, G{E-p) = 
^J^^dET{E) (-M|_^f)^^ ^jjgjg j(^) j, the zero- 
bias Fermi distribution in the electrodes and T{E) is the trans- 
mission coefficient for electrons of energy E passing from one 
electrode to the other For synthetic molecules such as oUg- 
othiophenes ITll r '^ El, - En and therefore T{E) exhibits 
Breit-Wigner 01611 or Fano resonances [17] in the vicinity of 
El, and En. For conjugated symmetric molecules with ex- 
tended wave functions along the backbones, resonant values 
of T(£') can be of order of unity and therefore at low-enough 
temperatures G{Ep) ~ 2e^//i whenever Ep aligns with a res- 
onance. In principle, this high conductance channel could be 



observable in conjugated single-molecule wires of arbitrary 
length, provided inelastic scattering is negligible. In practice, 
for most recently-studied families of molecules the measured 
conductance is much less than 2e^//i because Ep lies in the 
HOMO-LUMO gap. This means a serious handicap from the 
point of view of molecular circuitry, because the conductance 
will decrease exponentially as a function of length, and from 
the point of view of molecular sensing, because the movement 
of resonances produce by other molecules or environmental 
changes will have less effect in the middle of the gap. One 
approach to overcoming this problem is to introduce a third 
gate electrode, which shifts Ei,, En relative to Ep. In what 
follows we explore an alternative approach based on varying 
the electrode composition to yield a Ep close to En or Ei, 
and consequently high-conductance molecular wires without 
the need for a gate electrode. 

The idea of using alkaline electrodes to tune the work 
function of electronic contacts is well known in the organic 
LEDs (OLEDs) community [18, 19, 20]. In this paper, for 
the first time, we examine the effect of using alkaline elec- 
trodes for single-molecule electronics. We present a series 
of ab-initio calculations based on density functional theory 
(DFT) 12111 and non-equilibrium Green's functions (NEGF) 
1 2211 which determine the effect of tuning Ep. We exam- 
ine two molecules, the archetypal benzene- 1,4-dithiol (BDT) 
molecule that has been the subject of many experiments and 
theoretical works ||2|,|23|] (see Fig. ([T]i), and a longer molecule 
(namely l,4-bis[4- 24 (acetylsulphanyl)phenylethynyl]-2j6- 
dimethoxybenzene (molecule R3, shown in Fig. ([T])) 12411 
made of three aromatic rings, with the middle ring slightly 
twisted by the effect of oxygen-related side groups. We also 
examine molecules R5, R7 and R9 with 5, 7 and 9 aromatic 
rings, obtained by doubling, tripling and quadrupling the num- 
ber of phenyl rings in molecule R3, respectively. 



THEORETICAL METHOD 

We use the Smeagol code, ||9t| which interfaces NEGF to 
the SIESTA code 12611 and obtains self-consistently the den- 
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FIG. 1 : Real space projection of the density of states on an interval of 
0.2 eV around the Fermi level of the R3 molecule between (a) gold 
leads and (b) sodium leads. 



sity matrix and the transmission coefficients. We employ a 
basis set that includes polarization orbitals (SZP), which is 
good enough in this case |27]. To approximate the exchange 
and correlation we use the local density approximation (LDA) 
1 2811 . which works relatively well for organic molecules and 
noble or alkali metals. The real-space grid is defined by a 
plane-wave cutoff of 200 Ry. The leads are made of slices of 
fee (gold) or bcc (alkali) lattices grown along the (001) direc- 
tion, with 16 atoms per slice and periodic boundary conditions 
along the directions perpendicular to the transport axis, x and 
y. The system is also made periodic along z to avoid surface 
effects. The scattering region includes the molecule and five 
and six slices of the electrodes on the left and right parts of the 
unit cell, respectively. Each transport calculation has around 
250 atoms. 



PRELIMINARY CALCULATIONS 

To evaluate how the Fermi level aligns relative to the molec- 
ular levels we first calculate the energy states of the leads and 
each molecule separately. Since in SIESTA the energy ori- 
gin is arbitrary 129J it is necessary to use a common reference 
to compare the levels of different calculations. For this pur- 
pose we use a hydrogen molecule 10 A away from the slab or 
molecule, whose bonding orbital is taken as reference. The 
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FIG. 2: Relative position of the electrodes Fermi energies and the 
molecular HOMO and LUMO states calculated with LDA and taking 
as reference the bonding orbital of a hydrogen molecule away from 
the system. 



results are shown in Fig. (|2]l. As can be seen, all Fermi 
levels of the slabs sits inside the HL gap of the BDT. The 
Fermi level is smallest for gold and increases through the al- 
kali metals from Li to Cs. In the molecules, the HL gap is 
smaller in R3 than in BDT, as expected, and decreases as the 
number of rings increases towards R9. These results suggest 
that the conductance of molecular wires will change dramat- 
ically if gold is replaced by alkali leads. However, they de- 
scribe only isolated molecules and do not include bonding to 
the surface and charge transfer. After coupling to the slabs 
the molecular states will broaden into transmission resonances 
and the LUMO will shift upwards in energy towards E-p, due 
to charge transfer onto the molecule 11511 . To understand these 
effects we now perform complete transport calculations. 

Since in most of these systems it is not exactly known 
where the molecule would attach on the surface, it is neces- 
sary first to determine the most stable bonding configuration. 
We achieve this by calculating the energy as a function of the 
distance between a S-C6H4-SH molecule and a slab made of 
4 layers of fee Au or bcc Li, Na, K, R or Cs grown along the 
001 direction, as a function of the distance between the free 
sulphur atom and the surface. The results can be seen in Fig. 
(O. The same calculation was repeated for various positions, 
which include top, bridge and a 4-atoms hollow configura- 
tion. In the hollow configuration, steric interactions between 
some of the atoms of the molecule and the surface atoms on 
the corners of the square can give different energies. We also 
computed then the case where the molecule was rotated 45 
degrees along the axis perpendicular to the surface, but the 
results were practically the same. 

We found that in all cases the most stable configuration cor- 
responds to the hollow position. For short distances, however, 
the bridge can become the most favorable bonding configu- 
ration in the case of alkali leads. Its cohesive curve shows 
also two minima, as opposed to the one-minimum curves of 
other configurations. The second minimum appears when the 
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FIG. 3: Energy of a S-C6H4-SH molecule on top of a surface of fee 
gold (a) and bee lithium (b), sodium (c), potassium (d), rubidium (e) 
and cesium (f), as a function of the distance between the free sulphur 
and the surface for various configurations. 



TABLE I: Total charge gained by the molecule after connecting it to 
the electrodes obtained from the MuUiken populations. 



Au 



Li Na K Rb 



Cs 



A(3(e) 0.45 0.63 0.88 1.15 1.24 1.34 



sulphur atom crosses the line between both surface atoms and 
goes deep into the surface. We expect then that the bonding 
configuration undergoes a transition between the hollow and 
bridge states when the molecule is compressed. That would 
also affect the transport properties and could be detected in the 
experiments as a small peak in the conductance histograms. A 
general trend is that, as the alkali atomic number increases, the 
molecule can approach closer to the surface due to the larger 
lattice constant of the surface atoms. This, in turn, slightly 
increases the bonding distance between the sulphur and the 
surface atoms, something which is manifested by the broad- 
ening of the parabola in the cohesion curve from Li to Cs, 
and decreases the interaction between the molecule and the 
surface. 



ALKALI VS. GOLD ELECTRODES 

We first simulated the benzene molecule between gold elec- 
trodes and obtained the well known result that the Fermi level 
is close to the HOMO resonance 030113111 . When gold was re- 
placed by Li, however, the LUMO resonance moved towards 
the Fermi level and the conductance increased, as can be seen 
in Fig. (lUl. This was expected from the results of Fig. ^ 
due to the change in position of the Fermi level relative to 
the molecular orbitals. The same trend was also found for 
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FIG. 4: Transmission coefficients of a BDT molecule between fee 
gold (a) and bcc lithium (b), sodium (c), potassium (d), rubidium (e) 
and cesium (f). The dash-dotted line shows the difference between 
the transmission curves with and without molecule. 



the other alkali metals. However, for the heavier alkali ele- 
ments K, Rb and Cs there was direct transport between the 
electrodes even in the absence of the molecule. This effect 
was produced by the reduction of the distance between the 
leads due to the penetration of the molecule in the surface for 
large lattice constants and the increasing electronic delocal- 
ization of the outer electrons as the atomic number increases 
in these elements. This can be clearly seen in the lower panels 
of Fig. dUi where the direct transmission increases for high en- 
ergies, which coiTespond to less bounded electrons. In some 
cases the molecule acted also as an impurity reducing the total 
transmission between the leads. As a consequence, it was not 
easy to clearly distinguish the molecular contribution from the 
electrode contribution. For that reason, our attention moved to 
the R3 molecule, for which the electrodes are far enough. 

Fig. ^ shows the transport properties of the R3 molecule 
between gold and alkali leads. As can be seen, the transmis- 
sion coefficients are very similar in all cases and are charac- 
terized by the presence of two main resonances, correspond- 
ing to the HOMO and LUMO orbitals, and a large HL gap. 
In the case of gold, the Fermi level sits in the middle of the 
HL gap, which produces a very small zero-bias conductance. 
However, for the alkali elements the Fermi level moves to the 
LUMO resonance and dramatically increases the value of the 
conductance 1321 . The increase of the conductance can also 
be seen by plotting on real space the density of states (DOS) 
on an energy interval of 0.2 eV around the Fermi level, as 
shown in Fig. ([T]i. In the gold case the DOS of the leads does 
not penetrate into the molecule and electrons have to tunnel 
between both surfaces. When Ep moves towards the LUMO, 
DOS appears also in the molecule and makes it metallic. 

We found that the conductance increases from Li to K, but 
after that it decreases towards Cs. Such non-trivial behav- 
ior can be explained by taking into account the movement 
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FIG. 5: Transmission coefficients of R3 between gold electrodes (a) 
and alkali electrodes (b). 
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FIG. 6: Projected density of states on the first C atom, the S atom 
and the Au electrodes for a R3 molecule in the same configuration as 
in Fig. lO (a). 



of the Fermi level, charging effects and changes arising from 
bonding to different surfaces. As the alkali atomic number 
increases, the LUMO resonance moves downwards (see Fig. 
(|5]l) but as it starts to cross the Fermi energy, the charge on 
the molecule increases and the intralevel repulsion moves the 
states upwards again 1151 . The total amount of charge trans- 
ferred to the molecule increases with the atomic number of the 
alkali elements, as can be seen in Table (HJi, which is expected 
due to the increasing electropositivity and delocalization. At 
the same time, the lattice constant of the contact increases 
with the atomic number and therefore the space in the hol- 
low position becomes larger, which reduces the coupling be- 
tween the contact sulphur atom and the surface and decreases 
the width of the LUMO resonance. The first effect increases 
the conductance from Li to K, whereas the second decreases 
it towards Cs. 

We therefore expect a dramatic increase of the conductance 
if gold electrodes were replaced by alkali leads. Based on our 
calculations we predict that gold is one of the worst possi- 
ble choices for high-conductance molecular nanowires, since 
the conductance decreases exponentially as the length of the 
nanowire increases. 

From Fig. ^ it is possible to deduce some general proper- 
ties which can help to understand the influence of various fac- 
tors. For example it is noticeable that while the LUMO reso- 
nances in all cases have a similar width, the HOMO resonance 
is very sharp when the molecule is coupled to the alkali leads, 
as compared to the gold case, and its width decreases with the 
atomic number of the alkali element. This arises because the 
LUMO is delocalized along the molecule (see Fig. ([T]i) and 
it is better coupled to the contacts, whereas the HOMO is lo- 
calized in the middle of the molecule and is sensitive to any 
change in the coupling. Since the coupling to the alkali leads 
is smaller than the coupling to gold and decreases from Li to 
Cs, the width of the HOMO resonance is much smaller in the 
former and becomes sharper as the atomic number increases 



||3J. 



Additional information on the electronic structure and its 
influence ont the transport properties can be obtained by ana- 
lyzing the projected density of states (PDOS), which we show 
in Figs. ^ and (|7]i for the Au and Na configurations cor- 
responding to Fig. (|5]l. In both cases the transmission reso- 
nances associated with the HOMO and LUMO orbitals coin- 
cide with features in the PDOS of the C and S atoms related 
to the p_i orbitals, which are perpendicular to the transport 
direction and make on the molecular backbone the delocal- 
ized states near the Fermi level. Such orbitals penetrate into 
the molecule and generate the channels associated with the 
resonances. In the case of gold there is an additional pi con- 
tribution on the S near the HOMO, due to the strong inter- 
action with the surface atoms that move these orbitals up in 
energy, but it does not go inside the molecule because it al- 
most disappears in the C atom. The small transmission in the 
HL gap is produced by the absence of C states in this region, 
even though there are S states generated by hybridization with 
the surface states. The density of states of gold and sodium 
is finite in all the energy window and is mainly due to the s 
orbitals. Notice we plot the density of states per atom, so even 
in regions where it looks relatively small it has a finite contri- 
bution which allows the formation of a conductance channel, 
like for example in Fig. ^ between 0.5 and 1 eV. 



LENGTH DEPENDENCE 

The fact that the Fermi energy is close to the LUMO res- 
onance opens the possibility of having non-exponential de- 
pendence of the conductance as a function of the molecu- 
lar length, provided the Fermi level for longer molecules is 
still close enough. Such a possibility would allow the use of 
molecular wires with various lengths to connect different parts 
of a nanoscale circuit with almost no change in the conduc- 




FIG. 7: Projected density of states on the first C atom, the S atom 
and the Na electrodes for a R3 molecule in the same configuration as 
in Fig. ID (b). 



tance. We checked this possibility by using longer molecules 
with the same central unit and more benzene rings. To ensure 
that only the length was varied and no additional effects were 
introduced, the rings were made coplanar with the outer ones 
and the same number was added on both sides. The lengths of 
such molecules, which we call R5, R7 and R9, were 3.4, 4.7 
and 6.1 nm, respectively. 

It is well known that as the molecular length increases, the 
HL gap decreases due to the higher electronic delocalization 
along the molecule. The HL gap and the relative position of 
the HOMO and LUMO of all molecules can be seen in Fig. 
(|2]l. However, as the HL gap decreases, the value of the trans- 
mission coefficients inside the gap decreases due to the larger 
separation between the electrodes. This implies again that, if 
the Fermi level sits inside the gap, the conductance will de- 
crease exponentially as a function of the molecular distance. 
In Fig. ([8]l we show the transmission coefficients around the 
Fermi energy as a function of the distance for gold (a) and 
sodium electrodes (c). As can be seen in Figs, ©(a) and (b), 
which are plotted on a logarithmic scale, the transmission at 
the Fermi level decreases exponentially when the molecules 
are contacted to gold. Such exponential behavior has a beta 
of 0.16 A^^ and indeed proves that the quantum transport 
is in the tunnelling regime. As the molecular size increases, 
the LUMO resonance moves downwards and the distance be- 
tween levels decrease, but these effects are not enough to elim- 
inate the exponential behavior. 

However, in the case of sodium leads, shown in Figs. 
©(c), the Fermi level is located near the LUMO resonance 
for all molecules. This makes the length dependence non- 
exponential and almost linear, as can be seen in Fig. ©(d). 
Such evolution comes from an interplay between the down- 
wards shift of the LUMO resonance, charge transfer and 
changes in the coupling. The first molecule does not fit 
well into the linear trend due to the much bigger width of 
the LUMO resonance. Such an anomalous resonance is pro- 
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FIG. 8: Length dependence of the transmission coefficients and low 
bias conductance for molecules between gold electrodes ((a) and (b), 
respectively), and molecules between sodium electrodes ((c) and (d), 
respectively). 



duced by a very large interaction between the leads and the 
molecular central part, where this orbital has the biggest 
weight. At the same time, in order to keep the same charge 
on the molecule, this orbital is also shifted upwards. The 
LUMO resonances of longer molecules have a width which 
decreases and a position that increases slightly with the molec- 
ular length. The reduction of the width can be explained by 
taking into account that the weight of the LUMO decreases at 
the edges of the molecule and therefore the coupling between 
this orbital and the contacts tend to decrease as the length 
increases (the same, but more pronounced, happens to the 
HOMO orbital). The upwards movement of the LUMO res- 
onance looks counterintuitive, since as the molecular length 
grows the orbitals move downwards, but it can be explained 
by the upwards shift of the levels produced by the additional 
charge on the molecule. 



CONCLUSIONS 

The long-term goal of using high-conductance molecular 
wires as interconnects in sub- 10 nm electronics requires that 
HOMO or LUMO resonances coincide with the Fermi en- 
ergy of metallic electrodes so that the conductance does not 
decrease exponentially and all nanoscale electronic elements 
are connected with almost the same conductance. The pin- 
ning of the Fermi level at the LUMO is also important from 
the point of view of sensing, since the effect produced by 
other molecules or environmental changes, which is mani- 
fested mainly by movements of the resonances, will be much 
higher. In this paper we have shown that by analogy with cur- 
rent technology used in OLEDs, this can be achieved by using 
alkali metals [36] instead of the more commonly used gold. 
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